Nonuniversal dynamic conductance fluctuations in disordered systems 
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Sample-to-sample fluctuations of the time-dependent conductance of a system with static disorder 
have been studied by means of diagrammatic theory and microwave pulsed transmission measure- 
ments. The fluctuations of time-dependent conductance are not universal, i.e., depend on sample 
parameters, in contrast to the universal conductance fluctuations in the steady-state regime. The 
variance of normalized conductance, determined by the infinite-range intensity correlation Csit), is 
found to increase as a third power of delay time from an exciting pulse, t. Ci{t) grows larger than 
the long-range intensity correlation C2{t) after a time tq ~ y^(g)tD [to being the diffusion time, 
((j) being the average dimensionless conductance). 



I. INTRODUCTION 



The ability of a conductor to carry electric current is 
characterized by its conductance G equal to the ratio 
of the current carried through to the voltage across the 
conductor. G is finite because of scattering on impuri- 
ties and lattice defects, which prevents the electrons from 
propagating ballistically, along straight lines, transform- 
ing their trajectories into random walks with a step size 
£. The conductance G depends on the exact configu- 
ration of impurities inside the conductor. The ensem- 
ble average of the conductance, (G), can be expressed 
as (G) = {e^/h){g), where (...) represents the average 
over an ensemble of random impurity configurations, g is 
the dimensionless conductance, e is the electron charge, 
and h is the Planck constant. Deviations of conductance 
from its average value, 6G ~ G — (G), arise from the 
interference of electronic wavefunctions and, for this rea- 
son, are only important in the mesoscopic regime, i.e., 
when the electronic wavefunction is phase coherent. A 
prominent feature of electronic transport in mesoscopic 
systems is the universal conductance fluctuations (UCF): 
((5G^) ~ (e^/ft.)^, for any metal, independent of sample 
size or disorder strengtbii^iS, The UCF can be expressed 
in terms of g as {Sg"^) ~ 1. 

There is a remarkable analogy between electronic 
transport in mesoscopic systems and wave propagation 
in random media^i^i^. The dimensionless conductance g 
can be defined for classical waves as the transmittance 
T, that is, the sum over transmission coefficients Tab 
from an incident mode a to an outgoing mode 6, at a 
given frequency v, g = T = T^abTab- Here the modes a 
and b correspond to transverse eigenmodes of the empty 
waveguide and form a complete set of basis functions, 
sufficient to represent any wave field entering or leav- 
ing the waveguide. A given mode a is characterized by 
a transverse wave vector q^. The statement of UCF 
is then equivalent to "Eaba'b'iSTabSTa'b') ^ 1- The ori- 
gin of UCF can be traced to correlations of transmis- 
sion coefficients Tab of a random sample. Feng, Kane, 
Lee, and Stoned have expressed the intensity correlation 



function, Gaba'b' = {STabSTa'b'), as a series using l/{g) 
as the expansion parameter: Gaba'b' = Gi + G2 + G3. 
Gi ^ {TabY is the largest contribution to Gaba'b' but 
it is different from zero only when Aq^ = Aq;,, where 
^Qa = (la Qa' and similarly for Aqf,. It is therefore 
"short range" in both Aqa (at a given Aq^) and Aq?, (at 
a given Aq^). Its contribution to J^aba'b' {STab STa'b') is 
neghgible. G2 ~ {Tab)^/{g) is (5) times smaller, but it 
decays only algebraically with Aq^ and Aq^. For this 
reason, G2 is often referred to as a "long-range" corre- 
lation function. Its contribution to {ST"^) is still much 
less than unity. It is G3 ~ {Tab)'^ / {gY — the smallest of 
the three correlation functions and independent of a, a', 
b and b' — that explains UCF: var(T) ^ Y^aba'b'G^ ^ 1. 
Because G3 does not decay neither with Aq^, nor with 
Aqh, it is reffered to as an "infinite-range" correlation 
function. 

UCF have been extensively studied in mesoscopic sys- 
tems both experimentally^ and theoretically^ (see Refs. 
[l| andfsl for reviews). They have also been demonstrated 
in random media, with light^ and microwaves^, through 
the observation of the infinite-range (G3) intensity corre- 
lation. However, all studies of UCF yet conducted have 
only concerned the steady-state conductance. In this pa- 
per, we study fluctuations of the time-dependent con- 
ductance of a disordered system by means of diagram- 
matic theory and microwave pulsed transmission mea- 
surements. We show that the universality of conductance 
fluctuations is lost in the dynamic experiment: the fluctu- 
ations depend on the length of the sample and the mean 
free path The variance of normalized transmittance is 
observed to increase as a third power of delay time from 
an exciting pulse. A comparative analysis of the under- 
lying mesoscopic correlations allows us to identify a new 
characteristic time scale tq = \/ {9)^0, to being the dif- 
fusion time, at which G3{t) grows larger than C2{t). 
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II. THEORY 

In order to study the dynamic conductance fluctua- 
tions, we consider a pulsed wave propagating through 
a quasi-lD waveguide with static disorder, of cross sec- 
tion A and length L ^ £, parallel to the 2;-axis. This 
geometry also corresponds to the theoretical model of 
a disordered mesoscopic conductor. A pulsed excitation 
tpo{t) at a point po on the input surface Ain of the waveg- 
uide (z = 0) results in a complex intensity distribution 
/(p, t; Po) at the output surface ^out of the waveguide 
{z — L). To avoid confusion, we acknowledge the source 
position Po as an explicit argument of /. With a proper 
normalization of ^o{t)^ the transmittance T[t) is given 
by 



T{t)^ f <fp f d^poI{p,t-po), 



(1) 



and its variance is 



{dT{tf) - / d^p d^po 
X C{p,po;p',Po;t) 



d^p' j d^po' 

ut J A 



(2) 



where C(p, po] p' , p'o\ t) = {61{p, t; Po)SI{p', t; pf,)) is the 
time-dependent spatial correlation function of the inten- 
sity fluctuations, SI{p,t;po) = I{p,t;pa) - {I{p,t; po)) . 

In order to calculate C, we express the transmitted 
field ■ip{t) at a point r = {p, L} in terms of the exciting 
field ■ipQ{t') at ro = {po,0} with the help of the Green's 
function Q of Helmholtz equation: 



(3) 



V'(r,i;ro) = / dt'g{r,t;ro,t')Mt')- 

Using the definition of intensity, /(p, i;po) 
[■^(rjt; ro)P, we obtain 



Cip,po;p',p'o;t) = / l[dUK{r,ro,r',r'o;t,{U}) 

X Mtl)ro{t2)Mt3)ro{t4), (4) 

where the kernel K is given by the connected part of a 
product of four Green's functions, averaged over disorder: 

K = (g(r,i;ro,ti)g*(r,t;ro,i2) 

- (e(r,t;ro,ti)e*(r,t;ro,t2)) 

X {gir',t;rlt3)g%r',t;r'o,U)). (5) 

It is important to note that K can be written in the 
Fourier space as a correlation function of t/'s at four 
different frequencies. As a consequence, our calculation 
does not reduce to that of Ref. [y where the correlation 
function of steady-state conductances at two different fre- 
quencies was studied. Among many diagrams contribut- 
ing to K we are interested in those which do not vanish 
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FIG. 1: Two of the four diagrams contributing to the kernel 
K of Eq. © . The solid and dashed lines represent the Green's 
function Q and its complex conjugate Q* , respectively. Dot- 
ted hnes denote scattering events. H is the Hikami bos^*^, 
which describes an exchange of partners between 4 Green's 
functions. 



when both |p — p'| and |po — Pol exceed Aq, the wave- 
length in the medium at the carrier frequency vo- For 
a waveguide with a large number of transverse channels, 
N ^ I, such diagrams give the leading-order contribu- 
tion to {ST{t)'^). Similar to the steady-state case-'- i^, the 
diagrams of interest contain two pairs of Green's func- 
tions forming two ladders and exchanging partners twice 
on the way from ro and rg' to r and r'. We have identi- 
fied four such diagrams, two of them are shown in Fig. [T] 
They are similar to those encountered in the steady-state 
caseii^ii^, but include time-dependent Green's functions. 

When the diagrams for K are calculated, they are sub- 
stituted into Eq. (g]), which is used to calculate {6T{t)^). 
At this point it is necessary to assume a specific shape 
of the pulsed excitation. Here we consider a Gaussian 
pulse, |?/7o(t)P ~ exp(— t^/tp), and study the variance of 
the normalized transmittance s(t) = T{t)/{T{t)), 



var[s(i)] = 



{my 



(6) 



for various pulse durations tjr^. 

For pulse duration tp significantly greater than the dif- 
fusion time, tp ^ t]j ~ {L + 2zo)'^ /tt'^D, where D is the 
diffusion coefficient, zo i is the extrapolation length, 
and for short time delays from the center of the excit- 
ing pulse, t ^ tp, we recover the steady-state result, 
var[s(i)] = 2/15(.g)^. A particularly simple result is ob- 
tained for short pulses, tp <ti tzj, and long time delays, 
t ^ tu, when the average transmittance takes the form 
{T{t)) oc cxp(— i/i/)). Only the diagrams of Fig. [1] con- 
tribute to var[s(i)] in the leading order in tp/to- Keeping 
only terms that do not decay exponentially with t/tu, we 
obtain 



var[s(i)] — 



5V2 tp 

TT^y/ntD 



(7) 



This result is independent of absorption, if the pulse 
duration tp is much less than the absorption time ta, 
tp ^ ta- This condition can be readily met in prac- 
tice, opening interesting perspectives for probing scat- 
tering parameters of random media without necessarily 
knowing absorption inside the medium. For long pulses, 
tp > ti), all the four above mentioned diagrams con- 
tribute to ((5r(t)^). At long times t tp, var[s(f)] then 



takes the form of a cubic polynomial with the same coefh- 
cient in front of as in Eq. ([7]) and coefRcients in front of 
lower powers of t that have to be calculated numerically 
and will be presented elsewhere. 

To derive Eq. ([7]), we explicitly made use of the ex- 
ponential decay of the average transmittance with time: 
{T{t)) cx exp(— for t 3> t_D. Meanwhile, it is known 
that (T(t)) exhibits deviations from a simple exponential 
at long times t > t„ = ^/ (q)tn (weak localization correc- 
tions, see Refs.[llKilIl3i).' However, these deviations are 
of order 1/(5) for (g) ^ 1 and, therefore, taking them 
into account would yield corrections only of order l/{g)^ 
or higher in Eq. ([7]). In the present paper, we restrict 
ourselves to the calculation of var[s(t)] in the leading or- 
der in <C 1 and hence drop all higher-order terms. 
In addition, var[s(t)] = {5T {tf) / {T {t))'^ and corrections 
of similar type exist for the numerator {ST{t)'^) as well 
as for the denominator (r(i))^ of this ratio, leading to a 
non-exponential decay of both of them with time. As a 
result, a partial cancelation of these corrections might be 
expected in the final result. 

As is evident from Eq. ([7]), the fluctuations of the 
time-dependent conductance are not universal, in con- 
trast to the steady-state case. They explicitly depend 
on the mean free path, sample length, and other param- 
eters of the disordered sample. According to Eq. ([7]), 
at long times var[s(t)] can reach values well in excess of 
the steady-state result var(s) ~ 2/15 (g)^. Other conse- 
quences of Eq. ([7]) are the cubic growth of conductance 
fluctuations with time and their direct proportionality 
to the pulse duration tp, implying that T{t) is a self- 
averaging quantity in the limit of tp — > 0. 



III. EXPERIMENTAL SETUP AND SAMPLE 
CHARACTERIZATION 

Because the measurement of transmittance T requires 
the summation of Tab over all incoming and outgoing 
modes (or, equivalently, integration over input and out- 
put surfaces of the waveguide), it represents a real chal- 
lenge. In our microwave experiment, we make use of the 
relation between fluctuations of T and the intensity cor- 
relation function C3 discussed in the Introduction. Be- 
cause for vector waves, the notion of "mode" includes 
the polarization state as well, it is easy to show that 
if _L Ba' and _L e^/, the only contribution to the 
correlation function {STabSTa't') comes from C3 termS,; 
{STab{t)STa'b'{t)) = C'3(t). Here is the unit vector de- 
scribing the linear polarization state of the mode a. A 
link between var[s(t)] and C^it) is then readily found: 
var[s(i)] = Cz{t) / {Tab{t))'^ ■ Vector nature of electromag- 
netic waves leads to doubling of the number of transverse 
channels N as compared to the scalar case, and hence to 
doubling of {g) in Eq. 

Measurements of the correlation function of trans- 
mission coefficients {STab{t)6Ta'b'{t)) for _L Gq' and 
ef, _L eb' have been made for microwave radiation trans- 
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FIG. 2: (color online). Drawing of the experimental setup. 
Two microwave horns, S (source) and D (detector), are posi- 
tioned in front of and behind the random sample. The polar- 
izations of the incident and detected waves can be rotated by 
rotating the horns about their axes (z-axis). 



mitted through random mixtures of alumina spheres. 
Alumina spheres with diameter 0.64 cm and index of re- 
fraction 3.14 are embedded within Styrofoam shells to 
produce a sample with alumina volume fraction 0.03. 
The sample is contained within a copper tube with di- 
ameter of 4.4 cm and plastic end pieces. Linearly po- 
larized microwave radiation is launched and detected by 
conical horns placed 20 cm in front of and behind the 
sample, as shown in Fig. [51 Microwave field spectra 
are taken for cross-polarization orientations of both the 
source and the detector for 15 000 configurations pro- 
duced by briefly rotating the tube, in samples of length 
61, 76.2, and 91.4 cm. Measurements are made over the 
frequency range 18.7-19.7 GHz, away from sphere reso- 
nances and far from the localization threshold. At 19.2 
GHz, the number of transverse channels is TV = 54, with 
the effective refractive index of the sample, Ucs — 1.175. 
The transport mean free path estimated from Mie theory 
is ^ = 3.76 cm, giving {g) = ANe/3{L + 2zo) ~ 4.10, 3.33, 
and 2.81 for L = 61, 76.2, and 91.4 cm, respectively (we 
set zq = 2£/3). The temporal response to a Gaussian 
pulse of width tp peaked at t = is obtained by taking 
the Fourier transform of the field spectrum multiplied 
by a Gaussian envelope of width (2Trtp)~^ centered at 
1^0 — 19.2 GHz. The field of the time response is squared 
to give the transmitted intensity I{t). 

It is worthwhile to note that our experimental setup is 
similar to that of Ref. [3, except that the cross-polarized 
source and detector horns select the C3 correlation func- 
tion, whereas the measurements of Ref. [13 were domi- 
nated by Ci and C2. Indeed, one can define three dis- 
tinct transmission quantities: intensity transmission co- 
efficient Tab-, total transmission Ta = X^fj^ab trans- 
mittance T = X^ab ^afe- Fluctuations of Tab are described 
by Ci correlation function, whereas C2 and C3 govern 
fluctuations of Ta and T, respectivelyii^. Fluctuations of 
Tab and Ta were studied in Ref. [ij, whereas the focus of 
the present paper is on fluctuations of T. It is important 
to understand that because Ta is obtained from Tab by 
summing only over outgoing modes 6, whereas T includes 
the summation over incoming modes a as well, fluctua- 
tions of Ta and T are dominated by different scattering 
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FIG. 3: (color online). Time-dependent average transmission 
coefficient (in arbitrary units) of disordered microwave waveg- 
uides of three different lengths L = 61, 76.2 and 91.4 cm (solid 
lines). Polarizations of incident and detected modes a and b 
do not have any importance for this measurement because 
multiply scattered waves are completely depolarized. Dashed 
lines are theoretical fits to the data at short times (up to twice 
the peak arrival time) using diffusion approximation^^. Dif- 
fusion times to extracted from the fits are to = 17.5, 27.3 
and 39.6 ns, respectively. 



processes. 

To characterize our samples, we first estimate the ab- 
sorption time ta- To this end, we measure the decay 
rate of transmission in a sample of length L = 25.4 cm 
with copper end caps, weakly coupled to the measure- 
ment port (see Ref. |la for details of this measurement). 
We obtain ta = 39.7 ns. We then study the average trans- 
mission coefficient {Tab{t)). A fit of diffusion theory-^ to 
{Tab{t)) at short times (up to twice the peak arrival time 
of the pulse) yields diffusion times to — 17.5, 27.3, and 
39.6 ns for samples of lengths 61, 76.2, and 91.4 cm, re- 
spectively (see Fig. [3]). The fits in Fig. [3] exhibit notable 
deviations from the data at long times. This is expected 
from the previous workiiii^ii^ and is due to weak local- 
ization effects in our strongly disordered samples (see dis- 
cussion in Sec, mi . 
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FIG. 4: (color online). Variance of normalized transmit- 
tance, var[s(t)], in samples of three different lengths L for 
three different pulse durations tp (curves in the same order 
for the three panels). The diffusion time to extracted from 
the fit of diffusion theory to the average transmission coeffi- 
cient {Tab{t)) (see Fig. 131) is provided for each sample. Smooth 
lines plotted through the data points are the fitting curves of 
Eq. with {g) = 4.17(a), 3.21(b) and 1.85(c). Only the data 
for t > 3tD were used for the fitting. 



IV. RESULTS AND DISCUSSION 

The variance of normalized transmittance, var[s(i)], 
for the pulse durations tp = 0.9, 1.2, and 1.8 ns for each 
of the samples studied is shown in Fig. 01 The exper- 
imental data are fit with Eq. ^ using (g) as the only 
fitting parameter and imposing that (g) is the same for 
the three pulse durations for given L. The values of (g) 
found from the fit are 4.17, 3.21 and 1.85 for samples of 
lengths 61, 76.2, and 91.4 cm, respectively. These values 
differ slightly from the estimates, except for the longest 
sample. Going beyond the "two _ff-box" diagrams of Fig. 
[T][and hence, including higher-order terms in l/{g) in Eq. 
([7|] is necessary for a better quantitative agreement with 



the data for samples with small values of (g). 

To demonstrate that Eq. ([7]) correctly captures the key 
scaling properties of var[s(t)] for t ^ t]j, we introduce the 
dimensionless time r = t/to and the normalized pulse 
duration Tp — tp/tu, and plot var[s(t)]/rp versus r for 
L = 61 cm in Fig. [5] The three curves corresponding to 
the three different pulse durations fall on a single curve, 
demonstrating var[s(t)] oc Tp. The same data plotted on 
a log-log scale (inset of Fig. [5]) show that var[s(t)] cx t^. 
Thus, var[s(t)] cx TpT^, as predicted by Eq. ([7]). 

The power-law growth of var[s(i)] with time can be 
qualitatively understood from the following reasoning. 
In the path picture of wave propagation and according 
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FIG. 5: (color online). Scaling properties of var[s(t)] oc r^r"^, 
with Tp = tp/to and r = t/to, at long times, t 3> to, are 
demonstrated with the data for L = 61 cm and fp's as in Fig. 
|3ta). The theoretical (smooth) curve is for ((?) = 4.17. The 
inset shows the same plot on a log-log scale. 

to the diagrams of Fig. [U the value of var[s(t)] is pro- 
portional to the probability for two wave paths of length 
ct to cross twice inside the sample. The probability of 
a single crossing is p{t) ^ ct/NlJ^^. The probability 
that the two paths cross first during the time interval 
(0, t') and then a second time during the time interval 
(f, t) is p{t')p{t - t'). Var[s(i)] cx / {g)'^ is obtained by 
integrating over t' from to t. 

Keeping in mind the relation between var[s(t)] and 
C-i{t), we can use Eq. ([7]) to compare C:i{t) with 
the long-range intensity correlation function C2{t) oc 
i^/{9)){'tp/tD){t/tD), previously studied in Ref . [l3. Such 
a comparison suggests an important time scale, tg = 
V {9} to, intermediate between the Thouless time tTh ~ 
to and the Heisenberg time tn ^ {g)tD- This time scale 
appeared previously in the context of weak localization 
in classically chaotiot* and disordere d "'^^1^^ systems, but 



its role in the analysis of fluctuations of transport prop- 
erties has never been identified. Because for short pulses 
in the long-time limit we can write C2(t) ~ {tp/tq){t/tq) 
and C3(t) ^ {tp/tq){t/tqf' , the formally next-order con- 
tribution (in 1/(5) expansion) to the intensity correlation 
function, C^{t) becomes larger than C2(t) when t > tg. 
Even though this might indicate that some new physics 
come into play at such long times, var[s(t)] of Eq. ([7]) 
agrees well with the experimental data for t > tg, and 
even for t > tn, when the discrete mode structure of the 
disordered sample is expected to play a role^^. 



V. CONCLUSION 

In conclusion, we presented a theoretical and experi- 
mental study of fluctuations of the dynamic conductance, 
describing the response of a system with static disorder to 
a pulsed excitation. The variance of normalized conduc- 
tance is not universal and increases with delay time from 
an exciting pulse, t, reaching values well in excess of the 
steady-state variance. A comparative analysis of the un- 
derlying mesoscopic correlations C3{t) and C2{t) allowed 
us to identify a new characteristic time scale tg, being in- 
termediate between the Thouless and Heisenberg times, 
and to explain its role in dynamic statistics of transport. 
Our results apply to both electronic transport in meso- 
scopic systems and wave propagation in random media. 
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